Multiband Superconductivity in Spin Density Wave Metals 
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We study the emergence of multiband superconductivity with s- and d— wave symmetry on the 
background of spin density wave (SDW). We show that the SDW coherence factors renormalize the 
momentum dependence of the superconducting (SC) gap, yielding a SC state with an unconventional 
s-wave symmetry. Interband Cooper pair scattering stabilizes superconductivity in both symmetries. 
With increasing SDW order, the s-wave state is more strongly suppressed than the d-wave state. 
Our results are universally applicable to two-dimensional systems with a commensurate SDW. 

PACS numbers: 74.72.-h, 75.40.Gb, 74.20.Rp, 74.20.Fg 



Understanding the microscopic origin of thermody- 
namic phases with multiple order parameters (OP) is one 
of the central issues in condensed matter physics. This 
topic is particularly important for describing the com- 
plex phase diagram of many correlated metals such as 
the high-T c cuprates, ferropnictides, and heavy fermion 
compounds, in which it was suggested that superconduc- 
tivity coexists with a rich variety of other states, such 
as charge, spin, or orbital density wave states. Interest- 
ingly enough, in such coexistence phases, the presence 
of a density wave immediately leads to multiband super- 
conductivity due to the folding of the electronic bands. 
The relative phase of the superconducting (SC) order pa- 
rameters in multiband superconductivity was originally 
studied in Ref. and has attracted significant interest 
recently in the context of MgP>2 @ and the ferropnictides 
in which in-phase and out-of-phase locking, respec- 
tively, of the OPs has been found. 

The growing experimental evidence for the coexistence 
of superconductivity and a spin density wave (SDW) in 
the n-type (i.e., electron-doped) cuprates 0, ULEl, some 
heavy fermion [7[ and organic superconductors [8[ raises 
the question of how phase-locking occurs in such sys- 
tems. This question is of particular interest in the n-type 
cuprates due to a possible transition of the SC symmetry 
from (i-wave to s-wave with increasing doping. However, 
there remains experimental disagreement regarding this 
transition. While some measurements are consistent with 
a transition from a SC d x 2_ y 2-w&ve symmetry in the un- 
derdoped materials to either an s- or a (d-Hs)-wave sym- 
metry in the optimally and overdoped ones, other results 
suggest the existence of a d x 2_ y 2-w&ve OP with higher 
harmonics for the entire doping range [4|. 

In this Letter, we study the emergence of two-band 
superconductivity with d x 2_ y 2- or s-wave symmetry on 
the background of a commensurate SDW state with im- 
perfect nesting, using the n-type cuprates as a particu- 
lar example. We show that the SDW coherence factors 
renormalize the momentum dependence of the SC gap. 



This yields an s-wave OP which is unconventional, in that 
it acquires a 7r-phase shift between the two bands, line 
nodes along the boundary of the reduced Brillouin zone 
(RBZ), and changes sign between momenta connected 
by the SDW ordering momentum Q. In contrast, in the 
dj;2_ y 2-wave state, the OP is locked in-phase, with no 
additional line nodes. In both cases, superconductivity 
is stabilized by interband Cooper pair scattering. With 
increasing SDW OP, the unconventional s-wave state is 
suppressed more quickly than the d x 2_ y 2-wa,ve state. In 
the latter case, the vanishing of T c coincides approxi- 
mately with the disappearance of the first of the two 
Fermi surfaces (FS). To demonstrate the generality of 
our results for other materials, we consider both n- and 
p-type doping, leading to different evolutions of the FS 
in the SDW state. 

Our starting point for investigating the coexistence of 
superconductivity and SDW order is the Hamiltonian 

ko" k,k',cr 

+ C lk+qcr C p-q5 C pff C k(T (1) 

k.p.q.cr 

where (J (cko-) creates (annihilates) an electron with 
spin a and momentum k. We consider a two-dimensional 
system with the normal state tight-binding energy dis- 
persion £k = — 2< (cos k x + cos k y ) + At' cos k x cos k y — 
It" {cos2k x + cos2/c y ) — (i and hopping matrix elements 
t = 250 meV, t'/t = 0.4, and t"/t = 0.1. The chemi- 
cal potentials fi/t = —0.32 and fi/t = —1.09 describe the 
slightly underdoped n- and p-type cuprates, respectively, 
and the corresponding FS in the paramagnetic state are 
shown in Figs.[T](a) and (b). The second and third terms 
in Eq. ([TJ) give rise to a commensurate SDW and super- 
conductivity, respectively. While it is generally assumed 
that both phases emerge from the same underlying in- 
teraction, its renormalization due to vertex corrections 
gives rise to two different effective interactions, U and 
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FIG. 1: (color online) FS (green lines) for the n-type (a) and 
p-type (b) cuprates (x = 0.12) in the paramagnetic state. In 
the SDW state (Wo =0.1 eV), the a-band possesses electron- 
pockets (red) and the /3-band hole pockets (blue) . The RBZ is 
shown as a dashed line. The phase of the (a) unconventional s 
and (b) d x 2_ y 2 SC OPs are denoted by +/— . (c), (d) Angular 
dependence of the SC OP in the unconventional s-wave and 
d x 2 _ H 2-wave channels for n-type doping close to T c . 



V q , assumed to be independent of each other. 

In what follows, we study the coexistence phase in the 
limit T c <C Tsdwi and diagonalize the Hamiltonian in 
two steps. We first derive the electronic spectrum of 
the SDW state by decoupling the second term of Eq. (JTJ) 
via a mean-field (MF) approximation and then diago- 
nalize it together with the first term using an SDW Bo- 
goliubov transformation, introducing new quasi-particlc 
operators ak,/?k for the two resulting bands with dis- 

Here, 
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k , t- „„. W = 

U/2 cr( c k+Q o- c k,(r}sgncr is the SDW order parame- 
ter with ordering momentum Q = (%, tt), and e k = 
(ek ± e k+cj) /2. Due to imperfect nesting the a band ex- 
hibits electron pockets around (±7r, 0) and (0, ±7r) while 
the (3 band possesses hole-pockets around (±tt/2, ±7t/2), 
as shown in Figs. HJa) and (b) for the n- and p-type 
cuprates, respectively. The size of these pockets de- 
creases with increasing Wo- At Wo = W cr i ~ 0.19 eV 
(W cr i w 0.23 eV), the hole (electron) pockets disappear 
first for the n-type (p-type) cuprates, followed by the van- 
ishing of the electron (hole) pockets at W cr i > W cr i, as 
shown by the evolution of the density of states (DOS) 
in Figs. HJe) and (f). For W > W cr 2, the system is an 
antifcrromagnetic insulator (AFI). 

We next apply the SDW Bogoliubov transformation 
to the SC pairing interaction [third term in Eq.JT])], and 
subsequently performing a MF decoupling in the particle- 
particle channel, keeping only anomalous expectation 



values of the form 



- t t 
• a k,T a - 



k.| 



) and fj^lk i): an d their 



complex conjugates. The resulting MF Hamiltonian can 
be diagonalized by two independent Bogoliubov transfor- 
mations, yielding il k = 



k - y(K) +(AZr (7 = a, 13) as 



the energy dispersion of the two bands. 



The SC gaps, 



A k are determined self-consistently from two coupled 



gap equations given by (at T = K) 
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, The coupling of the OPs in the a and 

(3 bands results from momentum dependent interband 
Cooper pair scattering described by terms of the form 



tQ ^-k j.)<^k+q,T/3-k-q,i in the MF Hamiltonian. Be- 



low, we use V(k,k') = V s and V(k, k') = V^k^k'/^ 
with t/3k = cos k x — cos k y as the pairing interactions 
in the s- and d x 2_ y 2--w&ve channel, respectively. We 
checked that in the limit considered here, i.e, A sc <C Wo, 
feedback effects of superconductivity on the SDW order 
are negligible, and moreover, pairing terms of the form 
{a k +/3_ ]k ,) are irrelevant due to the FS mismatch be- 
tween the a and (3 bands. 

We begin by discussing the form of the SC OP. In the 
s-wave case, the SDW coherence factors entering the gap 
equation, 



Eq.([2]), can be factorized and one finds A k = 



with A£ = F c 



-F , and A£ follows via 

9 9\ A 7 , fw 
u l~ v l) T# tanh (# 



k ^k 
a <-» (3. Here, F 1 

with 7 = a, (3. Moreover, D k is unity if |-E k | < HuJd, 
and zero otherwise, with u>jj being the Debye frequency. 
The above form of the SC gap yields three important 
results. First, the s-wave gap is dressed by the SDW 
coherence factors, and hence acquires lines nodes along 
the RBZ boundary where u k — w k = 0. Second, there 
is a 7r-phase shift between the SC OP in the a- and j3- 
bands, as well as within the same band between momenta 
connected by Q, i.e., A k = — A k+ q. As a result, the 
s-wave symmetry of the SC OP is unconventional, as 
summarized in Figflja). Third, while Aq needs to be 
calculated self-consistently, its temperature dependence 
does not affect the momentum dependence of the SC OP. 



In the d-wave case, the SDW coherence factors can 
again be factorized, and one finds from Eq.([2]) A k = 
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FIG. 2: (color online) T c for unconventional s- and d-wave 
symmetries, as well as the DOS, as a function of Wo for (a), 
(c), (e) n- and (b), (d), (f) p-type cuprates. We set Wo = 
0.1 eV, and for the n-type (p-type) cuprates Va — 2 eV (Vd = 
1.2 eV) and V s = 0.6 eV (V s = 0.47 eV). 

<P]i(Aq + u k v k Aj) with = F Q a + Fg, Af = -A? = 

Ff - Ff. Here, = % £ p £>pV P £f tanh (§) and 

*7 = ^Ep^ZVp«p«P^tanh(g). Since |A?| < 

|A2|, the SC OP in the a and /3-bands are in phase 
with no additional line nodes, in contrast to the s-wave 
case. The resulting phase of the SC OP is shown in 
Figs. Gib). Note that the temperature evolution of the 
self-consistently determined Aq 1 can lead to changes in 
the momentum dependence of the d-wave OP. Finally, we 
find that for both SC symmetries, the SC OP evolves con- 
tinuously into that of the paramagnetic state, obtained 
in the limit Wq — > in which the folded parts of the FS 
disappear. In particular, for the s-wave case, the same 
sign of the SC gap is restored over the entire (large) FS. 

We next study the dependence of T c on the SDW or- 
der parameter, Wo- To this end, we linearize Eq.Q, and 
solve it iteratively in the RBZ on a 500 x 500 lattice, 
setting uj£> = 0.1 eV. This approach also yields the mo- 
mentum dependence of the SC gap at T c — + . A study 
of the temperature induced changes in the momentum 
dependence of the d-wave OP will be reserved for future 
work. In Figs. [51(a) - (d) we present T C (W )/T C (W = 0) 
as a function of Wo for the d x 2_ y 2- and s-wave channels 
and the n- and p-type cuprates. In both channels, T c 
decreases with increasing Wo- For the d-wave symmetry, 
we find that T c becomes exponentially suppressed once 
Wo exceeds W cr i, where the first FS pockets disappear 
in the pure SDW state. In order to understand this rapid 
decrease of T c around W cr i, we present in Fig. [3](a) the 
dependence of the effective intra- and interband interac- 
tion projected onto the s- or d-wave channel, defined via 

KTr = sL^p^p and v :?T r = eUO^p 

where ^ = 1 for the s-wave case, and 4> k — cos k x — cos k y 
for the d-wave case. For the d- wave case, V*V'j' ra decreases 
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FIG. 3: (color online) The projected inter- and intra-band 
effective interactions as a function of Wo- 



with increasing Wo, due to the vanishing of the SDW 
coherence factors in ££ < p for k — p = Q, in agreement 
with previous results [9(. The same argument, however, 
does not apply to V"^ L j er — 2 — V^fj ra which increases 
with increasing Wo, thus stabilizing the SC state. Once 
the FS of one of the bands disappears at W cr \ (indepen- 
dent of whether these are the electron or hole pockets), 
the channel for interband Cooper pair scattering, and 
hence T c , become rapidly suppressed and vanish slightly 
above W cr i due to the finite Debye frequency. While one 
might expect to find an exponentially suppressed, but 
non-zero T c as long as the system is metallic, i.e., for 
W cr i < Wo < W cr 2, its resolutions is currently beyond 
our numerical capabilities. This strong suppression of T c 
around W cr i provides an explanation for the experimen- 
tal observations that the emergence of hole pockets in 
n-type cuprates coincides approximately with the onset 
of superconductivity on the underdoped side [1, @, [l(| • 
For the s-wave case, T c decreases more rapidly than for 
the d-wave case, and becomes smaller than our numerical 
resolution already for Wo considerably smaller than W cr \. 
This behavior does not arise from a change in the DOS, 
which remains almost unchanged for Wo < W cr \ [see 
Figs.[2](e) and (f)], but is due to a decreasing magnitude 
of the effective interaction, V^rj ra — — V^j er , with in- 
creasing Wo , as shown in Fig. [31(b) . This decrease arises 
since the contributions to Eq.Q coming from the fold- 
ing of the original BZ onto the RBZ, are pairbreaking. 
In contrast, in the d-wave case, these contributions have 
the same sign as those coming from the RBZ, and thus 
are not pairbreaking. 

Finally, in Fig[T](c) and (d), we present the momentum 
dependence of the SC gaps for both symmetries along the 
electron and hole FS pockets as depicted in Fig. [T](a). 
In agreement with the above discussion, the s-wave OP 
is unconventional in that it possesses a 7r-phase shift be- 
tween the a- and /3-bands, and lines nodes along the RBZ 
boundary. In the cf^^^-wave case, the gap is nodeless on 
the a-band and possess a line node in the /3-band. While 
both the unconventional s- as well as the d-wave SC OP 
possess line nodes the slope with which the gaps increase 
as one moves away from the line nodes is considerably 
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FIG. 4: (color online) Schematic phase diagram of the n- 
type cuprates. The system is an antiferromagnetic metal with 
electron pockets for W cr i < Wo < W cr 2 (e-FS), and with 
electron and hole pockets for Wo < W cr \ (e+h-FS). 



larger in the s-wave case, where it is determined by the 
Fermi velocity, than in the <i-wave case, where it is set 
by the gap velocity [see Fig. [2(c) and (d)]. As a result, 
the DOS in both cases scales linearly at small energies, 
however, with a much smaller slope in the s-wave than in 
the d-wave case. In thermodynamic measurements, it is 
therefore difficult to distinguish between the unconven- 
tional s-wave OP discussed above, and a conventional, 
constant s-wave OP. We stress that the inclusion of the 
SDW coherence factors in the gap equation is crucial in 
determining the Wo-dependence of T c , and momentum 
dependence of the SC OP. Our results therefore lead to 
new properties of the coexistence phase, not discussed in 
earlier studies 11. 

The synopsis of our results, applied to the n-type 



cuprates, is shown in form of a schematic phase dia- 
gram in Fig. |U The system is an antiferromagnetic in- 
sulator for Wo > Wcri-, and becomes metallic with the 
emergence of electron pockets in the a-band at Wo = 
W cr i- While single-band superconductivity may occur 
for W cr i < Wq < Wcri-, the interaction necessary to 
achieve any appreciable T c is necessarily large rendering 
this possibility highly unlikely for both SC symmetries 
in real systems. However, the appearance of hole pock- 
ets in the /3-band at W cr i, allows for the emergence of 
two-band superconductivity with d-wave symmetry (blue 
curve). Upon further decreasing Wo, the current exper- 
imental results suggest one of two scenarios. First, su- 
perconductivity with an s-wave symmetry remains the 
subdominant OP, implying T^ < T c q even for Wo — > 
(red dashed curve). Second, if, as suggest by experi- 
ments, there is a transition of the SC symmetry from 
d- to s-wave (red solid curve), then this would natu- 
rally occur with decreasing Wq if for Wq — * one has 



TSo > T; 



C Q- The experimentally observed transition in the 



n-type cuprates between a state with well defined nodes 
to a state with large gapped phase space is consistent 
with this second scenario. Moreover, the observation of 
pair-breaking peaks in electronic Raman scattering which 
occur at similar energies in both Bi ff and P>2 g channels 
[l(J, would also suggest similar gap magnitudes, consis- 



tent with the form of the SC gap in the s-wave case (see 
Fig.[T](c)). Clearly, further, in particular, phase sensitive 
experiments arc required to clarify the SC symmetry over 
the entire doping range. 

In conclusion, we have studied the emergence of multi- 
band superconductivity with s- and d— wave symmetries 
in the presence of an SDW state. In both cases, the mo- 
mentum dependence of the SC OP is strongly renormal- 
ized by the SDW coherence factors. In the s-wave case, 
this leads to an unconventional OP which acquires line 
nodes along the boundary of the RBZ, and a ir phase shift 
between the a and (3 bands. For both symmetries, inter- 
band Cooper pair scattering stabilizes the superconduct- 
ing order. Moreover, with increasing Wo, the SC state 
with s-wave symmetry is more strongly suppressed than 
that with ci-wave symmetry. Our results are universally 
applicable to any (quasi-) two-dimensional system with 
commensurate SDW order in which the folding of the 
original BZ leads to separate FS pockets in both bands. 
Finally, we note that our results can be straightforwardly 
generalized to the coexistence of superconductivity with 
a commensurate charge or orbital density wave state [l3| . 
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